Momentum polarization: an entanglement measure of topological spin and chiral central charge 
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Topologically ordered states are quantum states of matter with topological ground state degeneracy and quasi- 
particles carrying fractional quantum numbers and fractional statistics. The topological spin 9„ = Inha is an im- 
portant property of a topological quasi-particle, which is the Berry phase obtained in the adiabatic self-rotation 
of the quasi-particle by 2n. For chiral topological states with robust chiral edge states, another fundamental 
topological property is the edge state chiral central charge c. In this paper we propose a new approach to com- 
pute the topological spin and chiral central charge in lattice models by defining a new quantity named as the 
momentum polarization. Momentum polarization is defined on the cylinder geometry as a universal subleading 
teiTn in the average value of a "partial translation operator". We show that the momentum polarization is a quan- 
tum entanglement property which can be computed from the reduced density matrix, and our analytic derivation 
based on edge conformal field theory shows that the momentum polarization measures the combination h^- 
of topological spin and central charge. Numerical results are obtained for two example systems, the non-Abelian 
phase of the honeycomb lattice Kitaev model, and the v = 1/2 Laughlin state of a fractional Chem insulator 
described by a variational Monte Carlo wavefunction. The numerical results verifies the analytic formula with 
high accuracy, and further suggests that this result remains robust even when the edge states cannot be described 
by a conformal field theory. Our result provides a new efficient approach to characterize and identify topological 
states of matter from finite size numerics. 



I. INTRODUCTION 

Since the discovery of integer and fractional quantum Hall 
states in 1980's''^, the research on topological states of matter 
has attracted tremendous experimental and theoretical inter- 
est. The recent discovery of topological insulators and topo- 
logical superconductors have greatly expanded the knowl- 
edge of topological statesj2. 1 In general, topological states are 
gapped zero-temperature states of matter which cannot be de- 
formed to trivial states without experiencing a quantum phase 
transition. Here the "trivial states" are reference states which 
have no quantum entanglement between any two different re- 
gions of space. Topological states of matter are characterized 
by topological properties such as gapless edge/surface states, 
ground state degeneracy, bulk quasi-particle excitations with 
fractional charge and fractional statistics. The fractionaliza- 
tion of quantum numbers and statistics can only occur in topo- 
logical states with ground state degeneracy, which are called 
topologically ordered states^'^. 

In a topologically ordered state, a given configuration of 
topological quasi-particles correspond to a finite dimensional 
Hilbert space. Braiding quasi-paiticles around each other 
leads to a unitary transformation in the Hilbert space, which 
is named as the fractional statistics of the particles^. When 
the braiding of different particles leading to non-commuting 
transformations, the statistics of the particle is non-Abelian, 
and otherwise it is Abelian. Besides the statistics, each quasi- 
particle a also has a fractionalized spin 6a which is a phase 
factor obtained by the state of the system during the self- 
rotation of a by 2n. The spins of topological quasi-particles 
are related to their statistics, because when two particles wind 
around each other by 2n and return to their positions, from 
far away, other particles cannot determine whether there are 
two particles braiding or there is only one particle (obtained 
by the "fusion" of the two particles) spinning around itself."^ 



The spin-statistics relation is a generalization of that for ordi- 
nary particles-/, e., fermions have half-odd-integer spin while 
bosons have integer spin. Therefore the spin values of quasi- 
particles are an important set of information that distinguishes 
different topological states. The calculation of topological 
spin is generically difficult except for algebraically defined 
topological field theories^ and some ideal model statesiSJ^. 
For microscopic lattice models with a topological phase, the 
topological spin of quasi-particles has only been computed 
numerically for several special cases: 1) The honeycomb lat- 
tice Kitaev model studied by calculating the braiding of quasi- 
particles'-'''*; 2) The fractional Chern insulators studied by 
modular transformations on a torus '^"'^. Such numerical re- 
sults are usually restricted by small system size and specific 
knowledge about the ground state wavefunction. 

For chiral topological states with robust chiral edge states, 
another fundamental topological quantity is the edge state chi- 
ral central charge. If the edge state is described by a chiral 
conformal field theory (CFT), it has a chiral central charge 
c, which determines the heat current Ie - jcT^ at a given 
temperature The central charge also appears in the gravi- 
tational anomaly of the edge'^ if the system is coupled to a 
gravitational field. Since the edge states of a topological state 
is only a well-defined one-dimensional state in energy below 
the bulk gap, it can at most be described by a CFT in the long 
wavelength limit. (As we will discuss later, there are exam- 
ples where the edge state is not Lorentz invariant even in the 
long wavelength limit, so that it's not a CFT at all.) 

In this paper, we propose a new numerical method for the 
calculation of topological spins and chiral central charge in 
candidate systems of topological states, which is significantly 
simpler than previous numerical methods, and can be gen- 
erally applied to identify topological states obtained in nu- 
merical approaches such as exact diagonalization and Monte 
Carlo. This method applies to a topological state defined on 
cylinder Since a cylinder is topologically equivalent to a 
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sphere with two punctures, each topological quasiparticle a 
corresponds to a ground state on the cylinder, which has quasi- 
particle a in each puncture. This is illustrated in Fig. [T](a). A 
rotation of the cylinder along the periodic direction y is equiv- 
alent to rotating the two particles simultaneously. By mapping 
to sphere we see that with reference to the direction normal to 
the sphere, one particle is rotated clockwisely while the other 
is rotated counter-clockwisely. Consequently, when the par- 
ticles have spin 9a, a 2n rotation of the cylinder should lead 
to a Berry phase of e'^" from one end and e"'*° from the other 
end, which cancels each other as expected, since this global 
rotation is a symmetry of the system. To avoid the cancela- 
tion, we need to define a "twist" of the cylinder where only 
the left half of the cylinder is rotated, as is illustrated in Fig. 
[T](b). For a lattice system, we do not have a way to twist the 
cylinder continuously, but have to twist discretely. If we con- 
sider the twist Fig. [T|(b) where half of the cylinder is twisted 
by one lattice constant, and there are Ly lattice sites in the y 
direction, such a twist should lead to a Berry's phase of e'^"/^' 
coming from the left edge of the cylinder with particle type a, 
without cancelation from the right edge. However, the twist 
looks like a violent operation since it induces a discontinu- 
ous jump at the boundary between left and right halves of the 
cylinder Naively this appeared to induce a non-universal cor- 
rection to the phase obtained during the twist, and thus make 
it impossible to measure the spin 0„. The key result of the cur- 
rent work is that the non-universal contribution is independent 
from topological sector a, which is shown in the following for- 
mula: 



Aa = {Ga\ \Ga) = CXp 



2ni 

—pa - aLy 



(1) 



with |Gfl) the ground state in the a sector, c is the central 
charge of the edge CFT, and a is a non-universal complex con- 
stant. The key property of this formula is that the first term 
is universal and contains the topological information ha and 
c, while the second non-universal term is independent from a. 
We name the coefficient pa of the first term as the momentum 
polarization, which has the physical meaning of the net y mo- 
mentum carried by the left edge in the ground state (in unit of 
Y^). We choose the unit which is the momentum quantum 

Ly Ly 

carried by a local bosonic excitation. Therefore Pa^'^ corre- 
sponds to a local excitation and the fractional part pa mod 1 is 
the topological contribution determined by the ground state. 
As we will show below, this result can be understood based 
on entanglement properties of topological states. 

The rest of the paper is organized as following. In Sec. 
we provide an analytic derivation of Eq. ([T) based on edge 
state conformal field theory and quantum entanglement prop- 
erties of the topological state. In Sec. |III] we verify this 
idea by numerical calculations in the honeycomb lattice Ki- 
taev model^**. In the non-Abelian phase of Kitaev model, we 
find that the spin of the non-Abelian any on cr v& Ba - f , in 
consistency with the expectation that this state is described 
by an Ising topological quantum field theory.^*^ In Sec. HVl 
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FIG. 1. (a) A global translation of the cylinder is equivalent to 
spinning two quasiparticles on a sphere to opposite directions, (b) 
The partial translation T';^ of the left part of the cylinder. In both 
panels, the red circles with arrows indicate the chiral edge states. 



this result is further verified in another topological state, the 
fractional Chem insulators——. The variational Monte Carlo 
wavefunction'^ of lattice Laughlin 1/2 state is confirmed to 
have central charge c - \ and quasiparticle charge Qa - f • 
Finally Sec. [V]is the conclusion and discussions. 



II. ANALYTIC RESULT 

In this section we will derive Eq. ^ based on the assump- 
tion that the edge state of the topological state is a chiral con- 
formal field theory. As will be seen later in Sec. |lll] the nu- 
merical results indicate that the applicability of formula ([T]) 
goes beyond the edge CFT assumption, but analytically the 
result can only be derived under this assumption. 

To begin with, we notice that the partial translation 
only acts on the left half of the system, so that we can trace 
over the right half and define the reduced density matrix 
PLa = trs (\Ga} <G„|) such that 



Aa={Ga\T^\Ga)^tlL{pLaT^) 



(2) 



The reduced density matrix pia is generically difficult to com- 
pute, but for topologically ordered states with chiral edge 
states, some long wavelength behavior of pi„ is known. In 
2008, H. Li and F. D. M. Haldane observed that the entan- 
glement spectrum (. i.e. the spectrum of the "entangle- 
ment Hamiltonian" defined as = -logp^ ) is qual- 
itatively the same as that of the edge Hamiltonian^^. This 
is proposed as a generic feature of the topological states, 
which has been verified numerically and analytically in var- 
ious physical systems such as general fractional quantum Hall 
wavefunctions^'''^^, free fermion topological states^**-^'' and 
Kitaev modeP*^, etc. A general analytic derivation of the 
entanglement-edge state correspondence has been proposed 
in Ref.^' for a generic topological state with chiral CFT edge 
states. Here we will briefly review the result of Ref^^ which 
will be applied to the current work. 
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The basic idea behind this derivation is that a cylinder can 
be considered as two cylinders A, B glued together, as is illus- 
trated in Fig. |2](a). Before the gluing procedure, each cylin- 
der has two edge states on the two boundaries, which in the 
long wavelength limit are described by a conformal field the- 
ory Hamiltonian Ha{B) - Ha{b)1 + HA(B)r, with I, r denoting the 
left and right edge of each cylinder A topological sector of the 
cylinder labeled by quasiparticle a corresponds to conformal 
blocks of the edge CFT formed by a primary field with scal- 
ing dimension hg and its descendants. The "gluing" is done 
by turning on a coupling between the two cylinders across the 
edge, and the coupling can be considered as a relevant cou- 
pling in the edge CFT. Using the property of CFT with rele- 
vant coupling, Ref.^' reaches the conclusion that the reduced 
density matrix of the edge state Ar is p^^l - Z^^' e"^'^"' ]^, 
which is a thermal density matrix of the chiral CFT restricted 
to the sector a. The left edge of A region is far from the B 
region, so that the left edge CFT stays in its ground state, with 
no entanglement with B region. This can be described by the 
same thermal density matrix with a temperature fii ^ oo. The 
total reduced density matrix of the left half cylinder (A region) 
is 

PLa = Pi7 ®pit = Z;' e-P'"--P'"-l , yS, = cx,, finite(3) 

as is illustrated in Fig. |2](b). Therefore Aa defined in Eq. dU is 
the average value of translation operator in a CFT with left 
edge at zero temperature, and right edge at finite temperature. 
In CFT the translation operator has the form 



(4) 



Here we have used the linear dispersion of the edge CFT 
Hi - vPi, Hr - -vPr with V the velocity of the CFT. Now 
we introduce the characters Xa(l) of the CFT^'^ which are 
defined as 



Xaiq) = IVaiq'^') = tra(q-"') = tl-„(^-«') 



(5) 



with tra the trace in the conformal block a. From Eq. (O and 
dUl, Aa can be expressed in the character by 



A, = 



X.(e-^"")xa(e-^''') 



(6) 



Now we consider a large but finite L,,, such that Vy6r ^ 
Ly <s; v/S/ = oo. The left edge is in zero temperature, so that 
the character is dominated by the ground state contribution 

^(q) = qh^-ij. Consequently = 



exp 



On the other hand, the right edge is in 

the high temperature limit since vj3r <K Ly. To compute the 
character in that limit we can make use of the modular trans- 
formation property of the character 



X. 



a{e-^''''-)^Y,^abXb(e'^n 



(7) 




FIG. 2. (a) A cylinder can be considered as two separate cylin- 
ders glued together. Each cylinder have chiral edge states (red and 
blue circles) and the gluing of the two cylinders couples the edge 
states at the interface, (b) After gluing the two parts of the cylinder 
are entangled through the interface, and the reduced density matrix 
of the left half cylinder is described by a thermal density matrix of 
the edge state CFT with finite "entanglement temperature" yS^' at the 
right edge, and zero entanglement temperature pj' = at the left 
edge (see text). 



with S at the modular S matrijs^. Taking e 
have 



Xa(e'''^"h^Yj^atXb[e-'-^) 



(8) 



Since the right-hand side is in low temperature limit, we can 
approximate I e 1 ~ e 24^ in the low temper- 

ature limit the trivial sector with hi =0 dominates all 
Other sectors with positive so that we hecve Xai^ ) - 
Saie~^. The same modular transformation can be applied 
to the numerator in Eq. (|6) and leads to Xa^^^^ ' ''^'j - 
•Sal 6^+^24 In summary we obtain 



exp 



-(ha--)- 
L, \ 24/ 



Ly- 



2nc 



(9) 



24 Vj3r [i + vPr) 

in the limit of Vy6,- <K Ly, which demonstrates Eq. ([TJ with 
a y^ (,-+1)3 ) non-universal, but independent from topolog- 



ical sector a. It should be noted that Rec : 



2nc 



>o. 



24 1)3,(1 +y2/32) 

which means Aa decays exponentially when the size Ly in- 
creases. 

For finite Ly/vf),-, the formula above obtains finite size 
corrections coming from two origins: i) the contribution of 
excited states to Xaiq) in low temperature limit. Xaiq) - 
q'^-^-k (1 + tiiaq + ...) with nia the number of states in the topo- 
logical sector a with momentum InjLy. ii) The contribution 
of other topological sectors b in the expansion ([8]l, which is 
proportional to q'''' . In the two corrections above, the most 
leading contribution comes from the ground state of the topo- 
logical sector with minimal scaling dimension. Denoting the 
minimal scaling dimension by we obtain 



Xa(f'^''\ 



Sa\e * 



l+O 



(10) 



and similarly for the numerator of Eq. (|6). Consequently the 
finite size correction to A„ is 



Aa = exp 



f (/.«-^)-aLyl(l.0(.-<^)^))(ll) 
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with a - 



2nc 



24 x'Prii + vPr) 



We see that the finite size correction to formula ([1} decays 
exponentially versus Ly/v/3r, which suggests that this method 
can be applied well to finite size numerical calculations. Such 
an exponential convergence behavior is confirmed in the nu- 
merical results presented in next section. 

Besides the finite size correction, other corrections to the 
formula ([T]i may occur in a realistic system due to deviations 
of the reduced density matrix p/,„ from the CFT behavior For 
example the edge state dispersion is not strictly linear, and the 
state counting may deviate from that of the pure CFT at high 
energy. Somewhat surprisingly, in the two physical systems 
we studied numerically, such nonuniversal corrections have 
not been found, and the result ([TJ remains valid even if the 
edge state deviates strongly from a CFT. 



III. NUMERICAL RESULTS ON HONEYCOMB-LATTICE 
KITAEV MODEL 

In this section we will calculate Ag numerically for the 
honeycomb-lattice Kitaev model"". Kitaev model is a special 
spin 1 /2 model which has a non-Abelian topological phase 
that we are interested in. For our purpose, we consider the fol- 
lowing Hamiltonian of the Kitaev model with nearest neigh- 
bor and second neighbor interactions: 

H ^ ^ JxO-ixO-jx + ^ JyCTiyO-jy + ^ JzO-izP-jz 
{ij)€x links {U)^y links {ij)^z links 

+Jnn ^ O-iyO-j.CTkx 



(12) 



(ijk)eA 



with X, y, z links the three types of nearest neighbor links on 
the honeycomb lattice shown in Fig. [3] (a). A in the J„„ term 
denotes the upper and lower triangles shown in Fig. |3](a), in 
which ijk runs counter-clockwisely around the triangle and j 
is at the upper or lower corner For simplicity, we only intro- 
duced one type of second neighbor coupling instead of includ- 
ing all possible terms obtained by 60 degree rotations. 

The unique property of the Kitaev model is that it can be 
solved by the following Majorana representation 



(13) 



in which a = x,y,z and 77, are four Majorana fermion op- 
erators defined at site /. The Hilbert space of 4 Majorana 
fermions is 4-dimensional. The two physical states of the spin 
at site / form a subspace of the Majorana fermion Hilbert space 
which is determined by the condition 

Di = WWiYim = -if^o^i^A = 1 (14) 
Using this representation the Hamiltonian is transformed to 

= - ^ ^ JaUijiriirij - Jnn ^ Uikiriirik 

a=x,y,z <(7>efllinks {ik)enn links 

with Uij = iylVj, for {ij) e a links 

Uik - UijUjk, for (ik) e nn links, (ijk) e A (15) 



Since [m^, Uki^ - for all dififerent links (ij) , (kl), the Hamil- 
tonian can be viewed as a free Majorana fermion Hamiltonian 
of 77, with Uij = +1 taking its eigenvalues. However the phys- 
ical states have to satisfy the constraint ( fT4l i. which can be 
obtained by a projection in the following form: 



Di + 1 



™«))>®|{«)> 



(16) 



with Di defined in Eq. ST4i . The link variable Uij can be 
viewed as a Z2 gauge field coupled to the Majorana fermion, 
and Di acts as a Z2 gauge transformation which transforms 
77, — > -77, and Uij -Uij for all links connected to /. The 
physical state is a superposition of all gauge equivalent 
configurations of gauge field Uij in direct product with the 
corresponding fermion state \^'f({u})}. Kitaev model has sev- 
eral phases including a gapless phase and Abelian and non- 
Abelian topologically ordered phases.^" In the following we 
mainly discuss the non-Abelian phase in which the 77, fermion 
has a band structure with a nontrivial Chern number C = 1 . 
There are three types of topological quasiparticles \ ,cr,^ in 
which 1 is the vacuum, ijj is the fermionic excitation of 77 Ma- 
jorana fermion, and cr is a Z2 flux of the gauge field M,y, with 
a Majorana zero mode of 77 Majorana fermion trapped in it. cr 
is the non-Abelian quasi-particle. Corresponding to the three 
quasi-particle sectors, there are three topological sectors for 
the system on a cylinder Depending on the flux of M,j gauge 
field around the cylinder, the fermion 77, has periodic or anti- 
periodic boundary conditions. As is shown in Fig. |3](b), the 
Majorana fermion has chiral edge states due to the topologi- 
cal band structure. When the Z2 flux in the cylinder is -1, the 
edge state has anti-periodic boundary condition and momen- 
tum eigenvalues k- 7^ (« + 5)- Denote the 77 fermion ground 

state in this sector as |G^), we have \G\) - P\G^f^ ® |{mo}) 
with P - n,(A+ 1)/2 the projection to gauge invariant states. 
Mq is a gauge configuration with flux - 1 in the cylinder For 
example the gauge choice in Fig. [3] (b) can be made, with 
Uij = - 1 for all red links across the dashed line, and Uij - 1 
otherwise. The 41 sector is the lowest energy excitation state 

\G^) - PflkfR-k I^f)®|{"o})' ^^^^fiRk the creation operator 
of the left and right edge states, respectively, and k - j-. The 
third sector cr corresponds to the periodic boundary condition 
\Ga) - P \G^f^ ® |{"o})' "0 1^ ^ guage configuration with -t-l 
flux in the cylinder, which can be taken as Uij - 1 for all {ij). 
It should be noticed that there are two degenerate states be- 
fore projection due to the Majorana zero modes, but only one 
survives in the projection. 

Due to the absence of quantum fluctuation of the Z2 gauge 
field, the entanglement entropy and entanglement spectrum of 
Kitaev model can be obtained rigorously^*'. Using the reduced 
density matrix obtained in Ref.^°, Ag in Eq. (|2]l can be re- 
duced to that of the free Majorana fermion: Aa - in^ppaTy^^. 
Here ppa is the reduced density matrix of the left half cylin- 
der, and tI:^ is the gauge covariant translation operator act- 
ing on the fermion 77,, which does a translation of the 77, 
fermions associated with a gauge transformation that trans- 
lates the gauge field configuration m^. This is like the mag- 
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FIG. 3. (a) Illustration of the definition of honeycomb lattice Ki- 
taev model in Eq. l ll2t . The x,y,z links are shown by the red solid 
lines, green dash lines and blue double lines, respectively. The sec- 
ond neighbor couplings are defined on the triangles shaded in orange 
color, with the labeling j, k of the triangle shown on the figure, (b) 
Illustration of the three ground states in the three topological sec- 
tors of the non-Abelian phase of Kitaev model. The right panels 
show the edge state energy spectrum with solid and hollow circles 
stand for occupied and unoccupied quasi-particle states. |Gi> and 
|G^^ coiTespond to the ground state and the lowest energy quasipar- 
ticle excitation state with fermion anti-periodic boundary conditions, 
while \Ga-) correspond to one of the degenerate ground states with 
fermion periodic boundary condition. The red dot in the right panel 
marks the Majorana zero mode, the occupation of which does not 
change the momentum polarization of IG,^). 



netic translation operators in a Landau level. More details 
on the reduced density matrix and the calculation of momen- 
tum polarization are given in Appendix |A] As a general re- 
sult for free fermions and free bosons^'*, pfa can be deter- 
mined by the two-point function (rjiri/j (due to the Wick the- 
orem), and the entanglement Hamiltonian has the quadratic 
form He = log Z„ - logpf„ = Y^ijeL iA/jnilj- The constant Z„ 
is determined by the normalization condition of the reduced 
density matrix. On the half cylinder, the entanglement Hamil- 
tonian can be diagonalized into the form 

HE-J]MfnJ''k-\\ (17) 

n,k ^ ' 

with k the y-direction momentum taking the values of j^n and 

Y (n + {n = 0,1,2, Ly- 1) for periodic and anti-periodic 
boundary conditions, respectively, fnk are the quasiparticle an- 
nihilation operators and ^„^. are the eigenvalues of the entan- 
glement Hamiltonian. In the momentum basis, the translation 
operator simply multiples a phase e''^ to each quasi-particle 
operator /„A , so that 



LF 



exp 



(18) 



Using Eq. (fTTl i and Eq. (fTSl l we obtain 

^ =z;'tr(e-'^^rf ) 



1 + e* 1 - 
+ 



^ -tanh — 

2 2 



(19) 



Here a — \,(t corresponds to the anti-periodic and periodic 
boundary conditions, respectively. Since the i/r state is an ex- 
cited state of the anti-periodic sector, it is not directly included 




FIG. 4. (a) The energy spectrum and (b) the entanglement spec- 
trum {^„i^ in Eq. dlVt ) of the Kitaev model, (c) The value of 
Ly (9„(Ly) + Q-L,.^ versus the system size L,. for a = 1 (blue x) and cr 
(red o) sectors. d„ = Imlog/i„. a is obtained from the slope of a lin- 
ear fitting of LyOaiLy) versus LJ. The dashed lines mark the theoreti- 
cal values given by Eq. (TJ Inlha - ^) for /i] = 0, /i,^ = -j^. (d) The 
deviation of 6a from the theoretical value 69a = 6a+aLy- 1^ [ha ^) 
versus Ly plotted in a log scale. The calculations are done for a sys- 
tem with L, = 40 along the x direction with zigzag edge and the 
parameters J;: = Jy = = I, J„n = 0.2. 



in the formula above. However, it is straightforward to show 
that = Aie'"^^' with the phase factor e™/^ contributed by 
the edge state fermion excitation f^. 

A honeycomb lattice on the cylinder can be defined with 
different orientations. Two different simple orientations are 
known as the "armchair edge" and "zigzag edge" in the 
graphene literature, as are illustrated in Fig. |3](a). We will fo- 
cus on the zigzag edge, which is convenient for the reason that 
will be clear in later part of this section. The numerical results 
of Act and Ai are shown in Fig. |4] To compare with the formula 
([T]l, we define OaiLy) - Im log Ag and do a linear fitting of LyOg 
versus LJ. Eq. ([T]) predicts Ly6a{Ly) - -l-maL^ + ljipa, so that 
the intercept of this fitting at Ly = gives the numerical value 
of ha and c. Fig. |4](a) shows 2npa - Ly6a(Ly) + IrnaL^, with 
the slope Ima obtained from linear fitting. The agreement 
with formula ([TJ is apparent. As is shown in Fig. |4](b), the 



-ImffLv + ^pa 



from the theoret- 



deviation 66a - OaiLy) 

ical value decays exponentially versus Ly, in agreement with 
the error estimation given by the CFT in Eq. (ITTt . 

An important question is how robust is the formula ([I}- 
Since the derivation given in Sec. is based on the CFT 
description of the edge states, it is not obvious whether the 
deviation of the edge state from CFT behavior will lead to 
connections to the universal l/Ly term in Eq. ([TJ. Interest- 
ingly, the numerical result above shows that the formula ap- 
pUes to the lattice model without any correction to the 1 /L,, 
term. To probe the robustness of the formula dl), we can mod- 
ify the Hamiltonian near the cut between left and right parts 
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FIG. 5. (a) The Hamiltonian of the Kitaev model is modified by 
turning off the second neighbor hoppings (puiple dotted lines with 
arrows) along the boundary between left and right parts of the cylin- 
der (purple dashed line), (b) The entanglement spectrum for this 
geometry, with quartic dispersing entanglement edge states, (c) The 
value of Ly [OaiLy) + q-L,,) and (d) the deviation of 6a from the theo- 
retical value for the modified geometry. The parameters are the same 
as in Fig. |4] 

of the cylinder. As is shown in Fig. |5] (a), we turn off the 
second neighbor hopping along the sites that are at the bound- 
ary between left and right parts. Such an edge term does not 
change the topology of the system, but changes the entangle- 
ment spectrum. As is shown in Fig. |5] (b), the edge state 
entanglement spectrum has a cubic dispersion Ek oc in the 
long wavelength limit, rather than linear. Therefore even in 
the long wavelength limit the entanglement Hamiltonian He 
is not a chiral conformal field theory. It is reasonable to be- 
lieve that the physical edge (left edge) is unaffected by such a 
change of dispersion, since it has zero entanglement temper- 
ature and stays in the ground state of the topological sector. 
However the entangled edge (right edge) is at finite entan- 
glement temperature, so that its contribution to momentum 
should be modified when the dispersion is changed. To our 
surprise, the numerical results shown in Fig. |5](c) (d) clearly 
shows that formula ([TJ still applies to this case. Such a robust- 
ness beyond edge CFT description indicates that the quantiza- 
tion of momentum polarization has a deeper origin from bulk 
topology. 



IV. PROJECTED WAVEFUNCTIONS FOR FRACTIONAL 
CHERN INSULATORS 

In this section we will calculate the topological spin 
and chiral central charge for a v = 1/2 lattice Laugh- 
lin stata^"— . These lattice Laughlin states, also known as 
the fractional Chem insulators, are obtainable from parton 
constructions^^-^"^, and we develop a Monte Carlo approach 
for calculating their momentum polarization. Based on Eq. 
([T]i, our numerical results for topological spin and chiral cen- 



tral charge are obtained for the v - Ijl lattice Laughlin state, 
which agree very well with the theoretical predictions from 
CFT. This suggests that momentum polarization is a very 
useful quantity for identifying chiral topological order from 
ground-state wave functions of interacting systems. 

The fractional Chern insulators are lattice analogies of frac- 
tional quantum Hall states^'"^'*. When bosons or fermions 
fractionally fill a nearly flat band with nonzero Chern number, 
incompressible quantum liquid states with fractional quan- 
tized Hall conductance can appear as stable ground states of 
a repulsive interaction Hamiltonian. In this work, we focus 
on a particular example of fractional Chem insulator, i.e. a 

V = 1 /2 lattice Laughlin state of hardcore bosons. Let us 
consider a N x N square lattice with single parton sites. To 
construct a v - 1/2 lattice Laughlin state, we use the par- 
ton construction in Ref. and split the hardcore boson at 
each site into two fermionic partons, bl - cj.|.c^. The physi- 
cal Hilbert space of hardcore bosons requires either two par- 
tons or no paiton at each site, coiTesponding to the pres- 
ence or absence of hardcore boson, respectively. To build 
a hardcore boson wave function from partons, a projector 
Pa is needed to remove those unphysical configurations with 
single-parton sites. Then, a hardcore boson state can be ob- 
tained by acting the projector Pc on the parton wave function 

It,,) = Pc I0«>T \<f>a)l- 

Now we assume that the fermionic partons are in the ground 
state of the following Chem insulator (See Fig. |6]l 

= X + ' X ^ikc]^Ck^ (20) 

Here the nearest-neighbor hopping integrals are equal to 
t along X direction and take the value f (-f) in odd (even) 
columns along y direction. The next-nearest-neighbor hop- 
ping integrals A,j are A (-A) if the hopping direction is along 
(against) the arrow in Fig. |6] In this work we focus on f = 1 
and A - 1/2. For periodic boundary conditions along both 
directions, the Hamiltonian ( |20] i has two bands with Chem 
number C = ±1. At half filling, the ground state of par- 
tons has a completely filled lower band. After projecting 
this parton band insulator state onto hardcore bosons, it has 
been shown"*" that the resulting projected wave function is a 

V = 1 /2 lattice Laughlin state with two-fold degenerate on a 
torus, and its topological sectors are related with the bound- 
ary conditions of partons: the identity sector (setnion sec- 
tor) is given by |Oi >= |0, > -i-|;r, >= |0,7r > -l-|7r,;r > 
(|(D, >= |0,0 > -\7T,0 >= |0,7r > -|7r,7r >), where |*P,,*P.v > 
is the projected wavefunction with periodic (^F^ - 0) or anti- 
period ("Pf = 7t) parton boundary conditions along x and sim- 
ilarly -Q,n along y. 

The parton construction of the v = 1/2 lattice Laughlin 
ground states on a cylinder has further complication due to its 
gapless chiral edge modes. To shed some light on the two- 
fold degenerate ground-state projected wavefunctions, imag- 
ine that we start on a toras and adiabatically lower all hop- 
ping amplitudes across the x boundary until they are much 
smaller than the chiral edge modes' finite size gap propor- 
tional to L As such a process involves no gap closing or 
level crossing, the wavefunction's topological properties, es- 
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FIG. 6. Illustration of the Chem insulator Hamiltonian l l20b . The 
two sublattices in the unit cell are labeled by A and B. The nearest- 
neighbor hopping integrals are / along x direction and odd columns 
along y direction, and -/ along even columns (dashed lines) along -y 
direction. The next-nearest-neighbor hopping integrals are A (-A) if 
the hopping direction is along (against) the arrow. 



pecially the eigenvalue of Wilson loops along y identifying 
the topological sectors should not change. Consider the spe- 
cial case when the chiral edge modes on the two edges of 
the cylinder have one exactly zero energy state each: Co-l and 
Co-R- On comparison with the open x boundary condition case, 
the small inter-edge hopping amplitudes will not modify the 
chiral edge modes, except for the coupling of the two zero 
energy states ci and cr in the form of 6{Yjiyc\_^^Ca-R +h.c}^. 
Here 5 > ( 5 < 0) is remnant from periodic (anti-periodic) 
boundary condition along x, respectively. In this way, the de- 
generacy between ci and cr is lifted, and the corresponding 
projected wavefunctions by filling all negative energy par- 
ton states are |0, - Pc (c^| - c^-f) - c^l) I'J?') and 
\n,^y') - Pc{c\^ h-c^JiO'), respectively, where 

|0') - ne„<o c|(q')c|(q;)|0) is obtained with all other parton 
states deeper in the valence bands filled. As previous knowl- 
edge about the torus geometry case would imply, the two topo- 
logical sectors on a cylinder correspond the symmetrization 
and anti-symmetrization of the 6 > Q and <5 < cases, which 
give: 

|0, > = Pa (4,4, + |0'> = Pc^T^i \'^') (21) 

|0, > = Pa (4,4, + 4,4,) lO') = Pacini, Id.') (22) 

In the last step of both equations we have used the particle 
hole symmetry of the Hamiltonian (l20l l. \<S>a) , a - \,s can be 
written in the equivalent form of lOa) = Pq \4'a)'\ \<l>a)i, with 
the two parton wavefunctions 10,,), , defined by 



e„<0 

|0,)^ = fk^^]^4(a)|O) 



V2 



(23) 
(24) 



e„<0 
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FIG. 7. Upper: The LyBc data of various values for the identity 
sector (blue crosses) and the quasiparticle sector (red dots) of a frac- 
tional Chem insulator. The intercept of the linear fits give the value 
of 2n(c/24 - ha). For the identity sector, the data is well consis- 
tent with linear fit (blue solid line) and suggests c = 1.078 ± 0.091; 
lower: the difference in Ly9„ corresponds to the measured values of 
h, = {LyOi - Lye,,)/2ji = 0.2521 ± 0.0063, the dotted line is the theo- 
retical value hs = 1 /4. 



The identification of topological sectors above is further 
confirmed by using a Monte Carlo algorithm'^ to calculate 
the overlap between a series of projected wavefunctions of the 

form \<S>(p,g)) = Pa \<f>(p,q))^ \4>{p,q\^ with 

k(P,?))^ = {pel + qcl^) n ^^(^)I0> (25) 

e„<0 

and p,q complex coefficients. The over- 

lap of the states with trial values (p,q) = 
(0,1), (1,0), (1/V2,+1/^), (1/V2,//V2) confirms 
that the only two linearly independent basis states are 
the two topological ground states {p,q) - (1,0) and 
(p,q)^(l/ V2, // V2), corresponding to Eq. (I23) and (l24t . 

To develop a Monte Carlo method to compute the momen- 
tum polarization, we note that the expectation value of the 
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twist operator for \(^a) can be express as 

{<i>a\T^\<^a) = V KH1>„)|'— r-f— (26) 

where |a) corresponds to a lattice configuration with 
hardcore bosons, and \T^a} a configuration with the posi- 
tions of hardcore bosons at the left half of the cylinder be- 
ing translated by one lattice spacing in y. Based on Eq. 
(|26] |. {'^a\Ty\^a) can be straightforwardly evaluated by treat- 
ing KfflOa)!^ as the probability for Monte Carlo updates and 
{<S>a\T^a}/{<i>a\a} as the Monte Carlo measurable. For the 
present v - 1/2 Laughlin state wavefunctions, these quanti- 
ties can be expressed in terms of determinant squares and up- 
dated efficiently with the inverse-matrix update techniques'*'. 
We emphasize that this Monte Carlo approach is also appli- 
cable for other chiral topological states from fermionic parton 
constructions'*^, whose wavefunctions take the form of deter- 
minants and/or Pfaffians. 

The results from this Monte Carlo method are shown in Fig. 
|7] According to Eq. [1] the intercept of the hnear fit of Lyfl^ 
versus for the identity sector corresponds to the measured 
central charge of the system c - 1.078 ± 0.091, whereas the 
difference in Ly0„ between the two sectors corresponds to Inh, 
so the measured topological spin is hs - 0.252 ± 0.006. These 
are well consistent with the known theoretical values for the 
model c - I and h = 1/4, and suggest Eq. ^ is also applica- 
ble to the parton constructed systems. 



V. DISCUSSION AND CONCLUSION 

In summary, we have proposed a new efficient method for 
computing the quasi-particle topological spin and chiral cen- 
tral charge for chiral topologically ordered states. By defin- 
ing the partial translation operator which translates the left 
part of the cylinder by one lattice constant, the combination 
of topological spin ha and chiral central charge c can be ex- 
tracted from a universal term ^ (jia - ^) in the phase of the 

ground state average values of in each topological sector 
h„ - ^ is the momentum polarization of the ground state in 
the topological sector a. We have derived the formula of cen- 
tral charge and topological spin (Eq. ([TJ) from the edge CFT 
form of the reduced density matrix, and also numerically ver- 
ified the formula in two distinct topological states, the honey- 
comb lattice Kitaev model and the v = 1/2 bosonic Laugh- 
lin state in fractional Chem insulators. Compared to previous 
numerical approaches which computes braid matrix of quasi- 
particles'^ '"* and modular S, T matrices'^"'^, the method we 
propose is more efficient. 

It is interesting to note that the momentum polarization de- 
fined on cylinder may be related to the characterization of 
one-dimensional symmetry protected topological states. If we 
view the cylinder as a one-dimensional (ID) chain, with the 
cylinder direction as an "extra-dimension", the system can 
be viewed as a one-dimensional system with a discrete 
symmetry. Here is the translation along y direction with 



N = Ly. In ID there is no intrinsic topological state but there 
are symmetry-protected topological states, which have non- 
trivial edge states at each end of the ID chain protected by 
global symmetries of the system. The type of topological state 
for a given symmetry group G is classified by the inequivalent 
projective representations of G, or equivalently by the second 
cohomology of the group H^(G, t/(l))^^. When the edge 
state carries a nontrivial projective representation of the group 
G, there will be a nontrivial series of group operations gi such 
that gNgN-[...gi = 1 but the sequential action of on the edge 
state leads to a nontrivial phase factor e'". This phase factor 
will be canceled the other edge so that the whole system is still 
invariant in G. Compared with the momentum polarization 
we see that the idea is very similar. A phase 9 - 2n(h - ^) 
is obtained by the left edge upon a 2?! rotation of the cylin- 
der, which is canceled by the right edge. However, the group 
cohomology of Z^ (and its continuum limit t/(l)) is trivial4^ 
Therefore if only the translation Z^ symmetry is concerned, 
the projective representation labeled by the fractional momen- 
tum h - can be linearized, and no nontrivial classification 
is obtained. More information is required besides translation 
symmetry to give a fundamental reason for the accurate quan- 
tization of ha, c in lattice models that is observed in the current 
work. Physically, the missing information is probably related 
to the locality of the model in 2D, which distinguishes a 2D 
model on the cylinder from a more generic ID model with Z^ 
symmetry. 

In the two example systems studied in this paper, the ground 
states in each topological sector is obtained from the knowl- 
edge to these states. However, it should be noted that in gen- 
eral our method does not require knowledge on topological 
sectors. If \Ga} , a = 1 , 2, n denote the n orthogonal ground 
states on the cylinder in an unknown basis, |Gq.) may be super- 
positions of different topological sectors |G„). An important 
observation is that the partial translation does not change 
the topological sector, so that {Ga\ \Gt} ^ 5 ah is diagonal. 
Physically, the reason is that topological sector a cannot be 
measured by local probes but can only be measured by cre- 
ating a pair of quasi-particles, taking one of them around the 
cylinder and then reannihilate them. We can choose the path 
of the quasi-particle so that it's far from the boundaries of 
the cylinder and the partition in the middle. As long as this 
is satisfied, such a quasi-particle motion is unaffected by the 
operation of partial translation T^. Therefore the topological 
sector has to remain the same after partial translation. Con- 
sequently, in a generic basis one just needs to calculate the 
matrix Kap - {Ga\ \Gi3^ and diagonalize it. The eigenval- 
ues are Aa with the form of Eq. ([T]l. 

Besides the examples given in the current work, the mo- 
mentum polarization can be applied a wide range of other 
topological states. The Monte Carlo method of computing 
the momentum polarization discussed in Sec. |IV] generally 
applies to other projected wavefunctions obtained by vari- 
ous parton constructions. In particular, this approach can 
be used to detect topological spin and chiral central charge 
in non-Abelian chiral topological states, which will be help- 
ful in identifying the topological order in more complicated 
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many-body wave functions. For instance, natural candidates 
for testing the momentum polarization scheme include SU(n)i: 
states'^ and SO(2n + l)i states^, whose trial wavefunctions 
and chiral edge CFTs are known. The momentum polarization 
results of these non-Abelian states will be reported in future 
works. 

Our approach may also have interesting applications in 
topological states described by density matrix renormaliza- 
tion group (DMRG)"*^. DMRG applies to all gapped ID states 
well, which is equivalent to finding a variational ground state 
with the form of matrix product states (MPS)^. MPS refers 
to states of the form |^> = 2(v,) [Ti'T^" ...T^^"]^^ |{v,)), with v; 

labels the physical states at each site, and TJ, a matrix with 
"internal" indices a,/3. The DMRG/MPS approach has been 
applied to 2D topological states on a cylinder^""^^, which can 
be viewed as a ID gapped state with a large number of states 
at each site, or one with a finite but long range interaction. 
Recently, the infinite DMRG, or infinite MPS, description has 
been studied for FQH states^^ and FCI states'^, which sim- 
plifies the variational state by assuming a translation invari- 
ant ansatz and taking the infinite length limit. Besides the 
MPS states obtained numerically, the model wavefunctions 
of FQH states obtained from CFT correlation functions have 
also been shown to have the infinite MPS form when defined 
on a cyUnderi^^i^ (with infinite internal dimension). The infi- 



nite MPS is particularly suitable for the study of the momen- 
tum polarization we propose, since it provides a description of 
the reduced density matrix using the internal auxiliary space 
(called the Schmidt states). In the long cylinder limit, the mo- 
mentum polarization can be directly computed from the ma- 
trices r'' at each site. An interesting open question is whether 
it is possible to obtain a generic proof of the universal quanti- 
zation of the momentum polarization hg- ^ from the infinite 
MPS approach, which, if exists, will be more rigorous than 
the CFT approach taken in the current work. Another inter- 
esting direction is to generalize the infinite MPS approach of 
Ref. |54 to cylinder FQH to FCI states by using the Wannier 
state representation^^"^?. 

We would like to note that Ref. [13 has studied the modular 
S , T matrices by a similar approach as Ref. Ts", and Ref. 113 
has independently studied the topological spin in an approach 
similar to ours, but with a continuous instead of discrete twist 
operator which is specially defined for cylinder FQH states. 
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Appendix A: Reduced density matrix and momentum polarization of Kitaev model 

We first review the structure of the reduced density matrix of Kitaev model obtained in Ref.-^". Without losing generality, in the 
following we assume L,, to be even. The first step is to separate the degrees of freedom to the left and right region. In the enlarged 
Hilbert space, the Majorana fermions 77, and most link variables u/j are akeady separated into the two regions L and R, and the 
only degrees of freedom that need special treating is the link variables on the links crossing the boundary. By pairing up the 
neighbor links across the boundary, we can redefine the link variables as is shown in Fig. |8]ln the 4 boundary sites 1,2,3,4, two 
link variables are defined as M12 - iTiT^, M34 = 'Tlyl- Now we define the new link variables - ij\Yy W42 - iylY2- It is easy 
to see that w,/s commute with each other and they both anticommute with m,/s. Also W13W42 - unu^A- In the 4-dimensional 
Hilbert space of the 4 Majorana fermions (before projection), choosing the eigenstates of M12, M34 we can find a basis in which 
M12 = cr^" i8> 1, M34 = 1 ® cr~, wi3 = o"-^ ® cr'', W42 - -cr^ ® cr^ ■ Therefore the eigenstates of wn, W42 are 

= ^ (ITi> + lit)) , = ^ (ITT) + IU» 
V2 V2 

I— ) = ^(ITT)-IU)), 1-+) = ^(ITi>-liT)) (Al) 
ITT) = (l++> + 1— » (A2) 

with the two ± labeling the eigenvalues of wn, W24 respectively, and T, i labels those of M12, M34. A state IjMiv}) with fixed values 
of Uij on all links can be written in the new basis. For simplicity we can always make a gauge choice so that = 1 for all links 
crossing the boundary. For such {M^j, we have 

hj^-^-'^"j:H^-j\),^h>>j\}, (A3) 

VV,;=±1 

A physical ground state has the form 

\G}=p(\GF{{uij}))®\{uij])) (A4) 
with P - Y\i ^^Y^- We write the 77, fermion ground state jCf ({mo})) in the Schmit decomposition form 

|g. ({«,})) = 2 ({<})) (A5) 

N 

with l*!*^^ and j*?* ) a set of orthogonal basis of left and right side rj fermions respectively, and q-ai the Schmit eigenvalues 
satisfying Yjn I'^'a'I^ = 1- Such a Schmit decomposition for quadratic fermion problem can be reduced to a simple single-particle 
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probleroas, the detail of which will be discussed later in this appendix. Thus the ground state is written as 



=±1 w 



(A6) 



in which we have written P - PlPr since the projection is local and separable to two subsystems. If we expand Pr - H/es ^r^' 
there are two terms which preserve the gauge field configurations ufj,Wij, which are 1 and YiieRDi- The second term does a 
gauge transformation to all sites t/,- — > -t/,- but preserves all the link variables. Therefore we can write 



K(K}))|{<'^-'7})> Z 



l + n,eRAUR 



(A7) 



and similar for P^. Here ^ denotes gauge equivalence. For two states Pr \^N\{^fjj)} |{"o' *^'v|/g with different boundary 
configurations Wjj, after the projection they are still orthogonal since it is not possible to do a gauge transformation to w,j 
without affecting Ujj on neighboring links. Also the number of orthogonal states in Eq. ( |A7| ) is the same for all configurations. 
Therefore the reduced density matrix is 

p, = Const. ■ J] J] \a,f P, [\% ({«(;.})) |{4, w,,}) J [{% {[uf^])\ {[u^, w,, 

w,j=±l N 

= Const. ■ ^ ^ la^l 



,2 i + n,ELA 



({4}))|{4>-;4K(K^)})|({4-^ 



l + n,ELA; 



{AS) 



The translation operator translates both t/,- and uij, wij. The global projection '^^2"^' only cares about the total fermion 
number parity, and commutes with T^. Thus 



A . Tr(r^p,) = Con.. ■ ^ ^ ^ (^^ 



^1n({4})) 



(A9) 



Here Ty'^ denotes the action of to 77 subspace, and Tu'fj is configuration obtained from the translation of u'fj by one lattice 
constant. For all gauge equivalent configurations the matrix element in Eq. iA9\ has the same value, since the gauge trans- 
formations to m' and Tu' cancels each other after being acted by T^. Therefore we can just take the representative configuration 



Uij and write 



A H Tr {T^pr) = Const. ■ ^ {[Tufj]} 



LWl^+UielD. 



,,=±1 N 



Now we do some more analysis to D,. We can write 



(MO) 



(All) 



i€L 



with the fermion number parity of 77 fermions. To write down this formula some suitable orientation convention should be 
made to the links since M,y - but the detail of this convention is not important here. The projector ''^n^'^'^ A jjj^posed the 
constraint YlieL F>i = 1 which requires - n<;/>eL Y\ij ^ij- Therefore in the sum over 2^>^^ possible configurations of w,j, 
half of them has the fermion projected to F^ - 1 and the other half projected to F,j - -1. Since the Schmit state j*?^ ({"0})) 
independent from the boundary variable w/j, the sum over the two projections simply give 



N 

Defining the gauge transformation V ({w^}) by 

^(K}44})K(K})) = K(K})) 



and the covariant translation operator 



Lrj 

y 



(A12) 



(A13) 



(A14) 
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we have 



withpif = J] \aM' ({«f,})> ({4})| 



(A15) 
(A16) 



To see the effect of T^^ more explicitly, we discuss the periodic and antiperiodic sectors separately. For the periodic sector 
with n<!7> M(7 = 1 around the cylinder, we can simply take the gauge choice u/j = 1 for all links. Since this gauge is already 
translation invariant, T^^ is the ordinary translation operator. The momentum eigenstates f„k has the wavefunction 



fnk - '^(f>nk(i.x)e 



with k - ^n, n e Z, and {ix, iy) the two-dimensional coordinate of / site. The translation operator acts as 

-LF _ ikj 



TLF 1 -rLF _ rpLF 1 f rpLt _ IK f 



(A17) 



(A18) 



The anti -periodic sector is a bit more complicated. We can choose a gauge with Ujj = - 1 on all links across a horizontal line, 
as is shown by the black solid line in Fig. [8](b), and = 1 on all other links. For convenience we label the horizontal chains 
parallel to the branchcut line by « = 1,2, L,, as is illustrated in Fig. |8](b) by the blue solid line. We choose the labeling such 
that the branchcut line is at the boundary, between n - 1 and n - Ly. In this gauge, the gauge transformation V ^jTM^ | , (wf^j) 
needs to translate the branchcut line by one lattice constant, which is defined by 



Vr]iV'' = -T]i, V / e chain 1 



(A19) 



Now if we define 



with yi e [1, Ly] and /: = |^ - ^) , n = 1, 2, Ly, the function e"'^^' has a jump between y,- - 1 and yi - Ly. Therefore 



Thus we see that the gauge transformation V correctly removes the additional minus sign and obtain 

rLF^l ^ 'rLF Jk r 
y Jnk^y - e Jnk 

Due to Eq. ( IA18I ) and ( IA22l l we see that T^^ in the particular gauge choice can be written in the form of Eq. ( fTsT i 



T';^ = exp 



n,k 



for both periodic and anti-periodic flux sectors. Using the fact that fj^i^fnk has eigenvalues 0, 1 we can write 



(A21) 



(A22) 



(A23) 



(A24) 



nk 



nk 



Using this equation and the entanglement Hamiltonian ( fTTI ). it is straightforward to obtain Eq. ([19). 

As the last part of this appendix we briefly review the free fermion reduced density matrix pif calculated in Ref.^'*. All 
multi-point functions (G| nf5i '7/ \G) of the free fermion ground state satisfies the Wick theorem, so that all correlation functions 
are determined by the two-point function 



Cij = {G\r]irij\G) - 6ij 



(A25) 
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(a) 




L ' R 



FIG. 8. (a) The definition of boundary variables Wjj on four sites 1,2,3,4. The boundary is drawn as the purple vertical dashed line, (b) 
Illustration of the gauge choice for anti-periodic boundary condition, with all ii/j = 1 except those crossing the branchcut line (thick black solid 
line). The sites on the lattice are organized into chains parallel to the branchcut line, shown by the blue thick solid lines. The chain labeling is 
given on the left. 



In particular, Cy for /, j e L in the left half cylinder determines all multi-point correlation functions in the left half cylinder. 
Therefore the reduced density matrix must i) also satisfy Wick theorem, and ii) reproduces the two-point functions. The first 
condition requires pif to have the quadratic (thermal) form 



Plf ^ exp 



i.jeL 



If we diagonalize He into the form given in Eq. (fTTT l 



the correlation function will be 



(A26) 



(A27) 



(A28) 



and corresponding Tr^pif /„j./j^), which should agree with the eigenvalues of C,j with the same eigenvectors. Therefore in the 
matrix level we have 



C = tanh/i =>/z= -[log(lH-C)-log(I-C)] 

This equation determines h and its spectrum from correlation function C. 

I 



(A29) 



